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Abstract 

A minimum cycle basis of a weighted undirected graph G is a ba- 
sis of the cycle space of G such that the total weight of the cycles in 
this basis is minimized. If G is a planar graph with non-negative edge 
weights, such a basis can be found in O(n^) time and space, where n 
is the size of G. We show that this is optimal if an explicit represen- 
tation of the basis is required. We then present an 0(n'^/^ logn) time 
and 0(n3/2) space algorithm that computes a minimum cycle basis 
implicitly. From this result, we obtain an output-sensitive algorithm 
that explicitly computes a minimum cycle basis in 0(n'^/^ logn -|- C) 
time and 0(n^/^ -|- C) space, where G is the total size (number of 
edges and vertices) of the cycles in the basis. These bounds reduce 
to 0(n^/^logn) and 0(n^/^), respectively, when G is unweighted. We 
get similar results for the all-pairs min cut problem since it is dual 
equivalent to the minimum cycle basis problem for planar graphs. 
We also obtain 0(n^/^logn) time and 0(n^/^) space algorithms for 
finding, respectively, the weight vector and a Gomory-Hu tree of G. 
The previous best time and space bound for these two problems was 
quadratic. From our Gomory-Hu tree algorithm, we obtain the fol- 
lowing result: with 0(n^/^logn) time and 0(n^/^) space for prepro- 
cessing, the weight of a min cut between any two given vertices of G 
can be reported in constant time. Previously, such an oracle required 
quadratic time and space for preprocessing. The oracle can also be 
extended to report the actual cut in time proportional to its size. 
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1 Introduction 



A cycle basis of a graph is a set of cycles that gives a compact representation 
of the set of all the cycles in the graph. Such a representation is not only 
of theoretical interest but has also found practical use in a number of fields. 
One of the earliest applications is in electrical circuit theory and dates back 
to the work of Kirchhoff [T7] in 1847. Knuth [18] used them in the analysis of 
algorithms. Cycle bases also play an important role in chemical and biological 
pathways, periodic scheduling, and graph drawing [15]. See also [U [H [6l [H 
[201 [231 [25]. 

In many of the above applications, it is desirable to have a cycle basis of 
minimum total length or, more generally, of minimum total weight if edges of 
the graph are assigned weights. The minimum cycle basis problem, formally 
defined below, is the problem of finding such a cycle basis. For a survey of 
applications and the history of this problem, see [H]. 

Let us define cycle bases and minimum cycle bases. Let G{V, E) be an 
undirected graph. To each simple cycle C in G, we associate a vector x 
indexed on where = 1 if e belongs to C and = otherwise. A set 
of simple cycles of G is said to be independent if their associated vectors are 
independent over OF {2). The vector space over this field generated by these 
vectors is the cycle space of G and a maximal independent set of simple cycles 
of G is called a cycle basis of G. Any cycle basis of G consists of m — + c 
cycles, where m is the number of edges, n the number of vertices, and c is 
the number of connected components of G [26] . 

Assume that the edges of G have real weights. Then a minimum cycle 
basis (MCB) of G is a cycle basis such that the sum of weights of edges of the 
cycles in this basis is minimized. The MCB problem (MCBP) is the problem 
of finding an MCB of G. 

The MCBP is NP-hard if negative weights are allowed [12]. The first 
polynomial time algorithm for graphs with non-negative edge weights was 
due to Horton [T^. His idea was to first compute a polynomial size set of 
cycles guaranteed to contain an MCB. In a subsequent step, such a basis 
is then extracted from this set using a greedy algorithm. Running time is 
O(m^n). This was improved in a sequence of papers [TJ [9l [21 [161 [HI [I] to 
0{'m'^), where u is the exponent of matrix multiplication. 

For planar graphs with non-negative edge weights, an 0(n^ log n) algo- 
rithm was presented in [12]. This was recently improved to 0(?i^) [I]. 

The quadratic time bound also holds for the following problem for pla- 
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nar graphs since it was shown to be dual equivalent to the MCBP for such 
graphs [12] (meaning that one problem can be transformed into the other 
in linear time): find a minimal collection of cuts such that for any pair of 
vertices s and t, this collection contains a minimum s-t cut. We refer to this 
problem as the all-pairs min cut problem (APMCP). 

We prove that quadratic running time for the two problems is optimal 
by presenting a family of graphs of arbitrarily large size for which the total 
length (number of edges) of all cycles in any MCB is 0(n^). 

We then present an algorithm with 0(n^/^ logn) running time and 0(n^/^) 
space requirement that computes an MCB of a planar graph implicitly. From 
this result, we get an output-sensitive algorithm with 0{n^^^ logn + C) time 
and 0(n'^/^ + C) space requirement, where C is the total size of cycles in 
the MCB that the algorithm returns. For unweighted planar graphs, these 
bounds simplify to 0{n^^'^logn) and 0(?t,^/^), respectively. Since the MCBP 
and the APMCP are dual equivalent for planar graphs, we get similar bounds 
for the latter problem. 

The weight vector of a weighted graph G is a vector containing the weights 
of cycles of an MCB in order of non-decreasing weight. Finding such a vector 
has applications in chemistry and biology [3]. From our implicit representa- 
tion of an MCB, we obtain an 0(n^/^ logn) time and 0(n^/^) space algorithm 
for finding the weight vector of a planar graph. The best previous bound was 
O(n^), obtained by applying the algorithm in [T]. 

A Gomory-Hu tree, introduced by Gomory and Hu in 1961 [TU], is a 
compact representation of minimum weight cuts between all pairs of vertices 
of a graph. Formally, a Gomory-Hu tree of a weighted connected graph G is 
a tree T with weighted edges spanning the vertices of G such that: 

1. for any pair of vertices s and t, the weight of the minimum s-t cut is 
the same in G and in T, and 

2. for each edge e in T, the weight of e equals the weight of the cut in 
G, defined by the sets of vertices corresponding to the two connected 
components in T \ {e}. 

Such a tree T is very useful for finding a minimum s-t cut in G since we 
only need to consider the cuts of G encoded by the edges on the simple path 
between s and t in T. Gomory-Hu trees have also been applied to solve the 
minimum fc-cut problem [23] . 
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For planar graphs, quadratic time and space is the best known bound for 
finding such a tree. The bound can easily be obtained with the algorithm 
in [1]. From our MCB algorithm, we obtain an algorithm that constructs a 
Gomory-Hu tree in only 0(?7,^/^ logn) time and 0{n^^'^) space. 

An important corollary of the latter result is that with 0{n^^'^ logn) time 
and 0{n^^'^) space for preprocessing, a query for the weight of a min cut 
(or max fiow) between two given vertices of a planar undirected graph with 
non-negative edge weights can be answered in constant time. Previously, 
quadratic preprocessing time and space was required to obtain such an oracle. 
The actual cut can be reported in time proportional to its size. 

The organization of the paper is as follows. In Section [2], we give some 
definitions and notation and state some basic results. We give the quadratic 
lower bound for an explicit representation of an MCB of a planar graph in 
Section [31 In Section |U we mention the greedy algorithm which has been 
applied in previous papers to find an MCB. Based on it, we present our 
algorithm in Section [5] and bound its time and space requirements. The 
corollaries of our result are presented in Section El In order for our ideas to 
work, we need shortest paths to be unique. We show how to ensure this in 
Section [71 Finally, we give some concluding remarks in Section [HI 

2 Definitions, Notation, and Basic Results 

In the following, G = (V, E) denotes an n-vertex plane, straight-line embed- 
ded, undirected graph. This embedding partitions the plane into maximal 
open connected sets which we refer to as the elementary faces (of G). Ex- 
actly one of the elementary faces is unbounded and we call it the external 
elementary face (of G). All other elementary faces are called internal. 

A Jordan curve partitions the plane into an open bounded set and an 
open unbounded set. We denote them by int{J) and ext{J'), respectively. 
We refer to the closure of these sets as int{J) and ext{J), respectively. 

We say that a pair of elementary faces of G are separated by a simple 
cycle C in G if one face is contained in int{G) and the other face is contained 
in ext{C). 

A set of simple cycles of G is called nested if, for any two distinct cycles 
G and G' in that set, either int{G) C int{G'), int{G') C int{G), or int{G) C 
ext{G'). A simple cycle G is said to cross another simple cycle G' if {G,G'} 
is not nested. 
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Figure 1: (a): A nested set B of five cycles Ci, C2, C3, C4, C5 defining five 
internal regions and an external region Roo{B) (white), (b): The region tree 
r{B) of B. 

For cycles C and C in a nested set B, we say that C is a child of C and 
C is the parent of C (w.r.t. B) if int{C) C int{C'). We also define ancestors 
and descendants in the obvious way. We can represent these relationships in 
a forest where each tree vertex corresponds to a cycle of B. 

For any cycle C G -B, we define internal region R{C, B) as the subset 
int{C) \ (Uj=i^,,,^fcmi(C.j)) of the plane, where Ci, . . . ,Ck are the children (if 
any) of C, see Figure [T](a). 

The external region R^oiB) is defined as the set \ (Uj=i^...^fcmi(Cj)), 
where Ci, . . . ,Ck are the cycles associated with roots of trees in the forest 
defined above. Collectively, we refer to the internal regions and the external 
region as regions. 

With Ci, . . . ,Ck defined as above for a region R (internal or external), 
we refer to the internal regions R{Ci, B) as the children of R and we call R 
the parent of these regions. Again, we can define ancestors and descendants 
in the obvious way. Note that the external region is the ancestor of all other 
regions. We can thus represent the relationships in a tree where each vertex 
corresponds to a region. We call it the region tree of B and denote it by 
T{B), see Figure [U^b). 

Note that for two cycles C and C in i3, C is a child of C if and only 
if R{C, B) is a child of R{C', B). Hence, the region tree T{B) also describes 
the parent/child relationships between cycles of B. 

The elementary faces of G belonging to a region R are the elementary 
faces of R. For each child Ci of i?, int{Ci) is called a non- elementary face 
of R. If R is an internal region R{C, B), the external face of R is the subset 
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ext{C) of the plane and we classify it as a non-elementary face of R. Col- 
lectively, we refer to the elementary and non-elementary faces of R as its 
faces. 

A cycle C in G is said to be isometric if for any two vertices u, v G C, 
there is a shortest path between u and v contained in C. A set of cycles is 
said to be isometric if all cycles in the set are isometric. 

The dual G* of G is the multigraph having a vertex for each elementary 
face of G and having an edge e* between two dual vertices for every edge e of 
G shared by the elementary faces corresponding to the two dual vertices. The 
weight of e* in G* is equal to the weight of e in G. We identify elementary 
faces of G with vertices of G* and since there is a one-to-one correspondence 
between edges of G and edges of G*, we identify an edge of G with the 
corresponding edge in G*. 

Assume in the following that G is connected. Given a vertex u ^ V, we 
let T{u) denote a shortest path tree in G with source u. The dual of T{u) 
is the subgraph of G* defined by the edges not in T{u). It is well-known 
that this subgraph is a spanning tree in G* and we denote it by T{u). The 
following lemma will prove useful. 

Lemma 1. Assume that for any two vertices in G, there is a unique shortest 
path between them in G. Let G he an isometric cycle in G and let u & V. If 
u E extiC) resp. u G intiC) then the elementary faces of G in int{C) resp. 
in ext{G) are spanned by a subtree ofT{u). If u E ext{C) fl int{G), i.e., 
u E C , then these two subtrees are obtained by removing the single edge of 
T{u) having one end vertex in int{G) and one end vertex in ext{C). 

Proof. Suppose that u G ext{G), see Figure [2j The subgraph of shortest 
path tree T{u) contained in int{G) is a forest. Since G is isometric and 
since shortest paths are unique, each tree in this forest contains exactly one 
vertex of G. This implies that the edges of G belonging to int{G) and not 
to this forest define a connected component in the dual of G. Since all these 
edges belong to T{u), it follows that the elementary faces of G in int{G) are 
spanned by a subtree of T{u), as desired. 

A similar argument shows that if m G int{C) then the elementary faces 
of G in ext{C) are spanned by a subtree of T{u). 

Finally, assume that u E G. There is at least one edge in T{u) with 
one end vertex in int{C) and one end vertex in ext{C) since otherwise, T{u) 
would be disconnected. There cannot be more than one such edge since that 
would contradict the first part of the lemma. This shows the second part. □ 
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Figure 2: If u G ext{C) then the subgraph of T{u) in int(C) is a tree. 

A Horton cycle of G is a cycle obtained by adding a single edge e to a 
shortest path tree in G rooted at some vertex r. We denote this cycle by 
C(r, e). For a subset V of V, we let HiV') denote the set of Horton cycles 
of G obtained from shortest path trees rooted at vertices of V. 

For any graph H, we let Vh and denote its vertex and edge set, 
respectively. If w : — M is a weight function on the edges of G, we say 
that a subgraph H of G has weight G M if J^eeEn '"^(^) — ^ ■ 

3 A Tight Lower Bound 

In this section, we show that there are planar graphs of arbitrarily large size 
for which the total length of cycles in any MCB is quadratic. This implies 
that the algorithm in [1] is optimal since it runs in O(n^) time. 

The instance Gn containing n vertices is defined as follows. Let wi, . . . , v„ 
be the vertices of G„. For i = l,...,n — 1, there is an edge Cj = (fi, Wj+i) of 
weight 0. For i = 1, . . . , n — 2, there is an edge e[ = (f i, fi+2) of weight 1. 

Since G„ has m = 2n—3 edges, any MCB of Gn consists of m—n+1 = n—2 
cycles. In such a basis, every cycle must contain at least one of the edges e^, 
i = 1, . . . ,n — 2. Hence, the cycles in any MCB of (?„ have total weight at 
least n — 2. 
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1. initialize S = 

2. for each simple cycle C of G in order of non- decreasing weight, 

3. if there is a pair of elementary faces of G separated by C and not by 
any cycle in B, 

4. add C toB 

5. output B 

Figure 3: The generic greedy algorithm to compute the GMCB of G. 

For i = 1, . . . , 77, — 2, let Cj be the cycle containing edges ei, . . . , Cj+i, in 
that order. It is easy to see that the set of these cycles is a cycle basis of G. 
Furthermore, their total weight is n — 2 so by the above, they must constitute 
an MCB of Gn- In fact, it is the unique MCB of Gn since in any other cycle 
basis, some cycle must contain at least two weight 1 edges, implying that the 
total weight is at least n — 1. 

The cycles in the unique MCB of Gn clearly have quadratic total length. 
This gives the following result. 

Theorem 1. There are instances of planar graphs of arbitrarily large size 
n for which the cycles in any MCB for such an instance have total length 

In Section we show how to break the quadratic time bound by com- 
puting an implicit rather than an explicit representation of an MCB. 

4 The Greedy Algorithm 

In the following, G = {V, E) denotes an ?7,-vertex plane, straight-line embed- 
ded, undirected graph with non-negative edge weights. We may assume that 
G is connected since otherwise, we can consider each connected component 
separately. We require that there is a unique shortest path in G between any 
two vertices. In Section [TJ we show how to avoid this restriction. 

The algorithm in Figure E] will find an MCB of G (see [H [H]). We call 
this algorithm the generic greedy algorithm and we call the MCB obtained 
this way a greedy MCB (GMCB) (of G). We assume that ties in the ordering 
in line 2 are resolved in some deterministic way so that we may refer to the 
cycle basis output in line 5 as the GMCB of G. The following two results are 
from [12]. 
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Lemma 2. The GMCB is isometric and nested and consists of Norton cycles. 

Lemma 3. For every pair of elementary faces of a plane undirected graph H 
with non-negative edge weights, the GMCB of H contains a minimum-weight 
cycle C in H that separates those two faces. Cycle C is the first such cycle 
considered when applying the generic greedy algorithm to H . 

Our algorithm is essentially the generic greedy algorithm except that we 
consider a smaller family of cycles in line 2. The main difficulty in giving an 
efficient implementation of the greedy algorithm is testing the condition in 
line 3. Describing how to do this constitutes the main part of the paper. 

5 Divide-and-Conquer Algorithm 

The family of cycles that we pick in line 2 of the generic greedy algorithm is 
obtained with the divide-and-conquer paradigm. 

To separate our problem, we apply the cycle separator theorem of Miller [22] 
to G. This gives in linear time a Jordan curve J intersecting 0{\/n) ver- 
tices and no edges of G such that the subgraph Gi of G in int{J) and the 
subgraph G2 of G in ext{J') each contain at most 2n/3 vertices. We let Vj 
denote the set of vertices on J' and refer to them as boundary vertices of G. 

As in G, we assume that shortest paths in Gi and G2 are unique. In 
Section [71 we show how to avoid this assumption. 

For i = 1, 2, let Bi be the GMCB of d. Let B'^ be the subset of cycles of 
Bi containing no vertices of Vj. 

Lemma 4. With the above definitions, B[ UB'2U T-C{Vj) contains the CMCB 
ofG. 

Proof. Let B be the GMCB of G and let C be a cycle of G not belonging to 
B[UB'^U n{Vj). We need to show that G ^ B. 

By Lemma [2], we may assume that G is isometric. Furthermore, we may 
assume that it does not belong to B1UB2 (since otherwise, it would belong to 
Si U ^2 \ {B[ U B2) and hence to T-C{Vj) since it is isometric and since shortest 
paths are unique). Since G ^ Ti-iVj), G does not contain any vertices of Vj 
so it belongs to Gi, where i G {1,2}. In particular, it is considered by the 
generic greedy algorithm in the construction of Bi. 

Since G ^ Bi, Lemma [3] implies that every pair of elementary faces (/i, 
of Gi, where fi C int{G) and /2 C ext{G), must be separated by some cycle 
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of Bi having smaller weight than C (or a cycle having the same weight as C 
but considered earlier in the generic greedy algorithm). We claim that this 
statement also holds when replacing Bi by B and Gi by G. If we can show 
this, it will imply that C is not added to B by the generic greedy algorithm. 

So let (/i,/2) be a pair of elementary faces of G with /i C int{C) and 
/2 C extiC). Either fi or /2 is an elementary face of Gi since either J C 
ext{G) oi J d int{C). Assume w.l.o.g. that /i is an elementary face of Gi. 

If /2 is also an elementary face of Gi belonging to the same connected 
component K of Gi as /i, the above implies that /i and /2 are separated by 
some cycle C G Bi having smaller weight than C . Since G' is also considered 
by the generic greedy algorithm when constructing B, it follows that /i and 
/2 are separated by a cycle in B having weight smaller than that of C, as 
desired. 

Conversely, if /2 is not an elementary face of Gi belonging to /T, /2 must 
be contained in the external elementary face fx of K. By Lemma [3], there is 
a cycle of Bi which is shorter than G and which separates /i and fx- This 
cycle also separates /i and /2 and it follows that /i and /2 are separated by 
a cycle in B having weight smaller that that of G. 

The above shows that G ^ B, completing the proof of the lemma. □ 

Lemma H] suggests the following divide-and-conquer algorithm for our 
problem: recursively compute GMCB's Bi and B2 of Gi and G2, compute 
niVj), and extract from B[ UB'^U HiVj) the GMCB of G by applying the 
generic greedy algorithm to this smaller set of cycles. Pseudocode of this 
algorithm is shown in Figure H] (it is assumed that a brute- force algorithm is 
applied to find the GMCB of G when G has constant size). We call it the 
recursive greedy algorithm. 

We will show how to implement the top-level of the recursion in Oin^/"^ logn) 
time and 0(?7,^/^) space. Since each step of the recursion partitions the graph 
into two subgraphs of (almost) the same size [22], it will follow that these 
bounds hold for the entire algorithm. 

Since the algorithm constructs the GMCB, B is isometric and nested at 
all times. Thus, B represents a set of regions that change during the course 
of the algorithm. More specifically, when the algorithm starts, B = ^ and 
there is only one region, namely the external region R^{B). Whenever a 
cycle G is added to B in line 5, the region R containing G is replaced by two 
new regions, one, contained in int{G) and one, R2, contained in ext{C). 
We say that G splits R into -Ri and R2. We call -Ri the internal region and 
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1. recursively compute GMCB's Bi and B2 of Gi and G2, respectively 

2. initialize S = 

3. for each cycle C E B[U 'HiVj) in order of non- decreasing weight, 

4. if there is a pair of elementary faces of G separated by G and not by 
any cycle in B, 

5. add GioB 

6. output B 

Figure 4: The recursive greedy algorithm to compute the GMCB of G. For 
i = 1, 2, B'^ is the set of cycles of Bi not containing any vertices of T-iiVj). 

R2 the external region (w.r.t. R and G). Figure [5] gives an illustration. 

The following lemma relates the test in line 4 to the two regions generated 
by the split. 

Lemma 5. The condition in line 4 in the recursive greedy algorithm is sat- 
isfied if and only if G splits a region into two each of which contains at least 
one elementary face. 

Proof. Let R be the region containing G and suppose that G splits R into 
Ri and R2. 

Consider two elementary faces of G separated by G. No cycle of B sepa- 
rates them if and only if the two faces belong to the same region. Hence, the 
condition in line 4 is satisfied if and only if G separates a pair of elementary 
faces both belonging to R. The latter is equivalent to the condition that 
there is an elementary face in Ri and an elementary face in i?2- CH 

5.1 Contracted and Pruned Dual Trees 

Lemma [5] shows that if we can keep track of the number of elementary faces 
of G in regions during the course of the algorithm, then testing the condition 
in line 4 is easy: it holds if and only if the number of elementary faces of 
G in each of the two regions obtained by inserting G is at least one. In the 
following, we introduce so called contracted dual trees and pruned dual trees 
that will help us keep track of the necessary information. First, we need the 
following lemma. 
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Figure 5: Adding a cycle C to B splits a region R into internal region Ri 
and external region R2. 

Lemma 6. Let H be a plane graph with non-negative edge weights and as- 
sume that shortest paths in H are unique. Let C he an isometric cycle in H 
and let P he a shortest path in H hetween vertices u and v. If hoth u and 
V helong to int{C) then P is contained in int{C). If hoth u and v helong to 
ext{C) then P is contained in ext{C). 

Proof. Suppose that u, v G int{C) and assume for the sake of contradiction 
that P is not contained in int{C). Then there is a subpath P' of P between 
a vertex u' E C and a vertex v' E C with all interior vertices belonging to 
extiC). Since C is isometric, there is a shortest path P" contained in C 
between u' and v' . But P' is also a shortest path between u' and v' . Since 
pi _^ pii ^ ^Yiis contradicts the uniqueness of shortest paths in H. 

A similar proof holds when u,v E ext{C). □ 

For a region R and a boundary vertex v belonging to i?, the contracted 
dual tree Tfi{v) is the tree obtained from dual tree T{v) by contracting each 
edge (n, u'), where u and u' are elementary faces in G both contained in the 
same non-elementary face of i?, see Figure [61 

An important observation is that there is a one-to-one correspondence 
between the vertices of Tr{v) and the faces of R. We assign the colour white 
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Figure 6: Contracted dual tree Tii{v) is obtained from T{v) by contracting 
edges between elementary faces belonging to the same non-elementary face 
(bold edges and white interior) of R. For this instance, applying the pruning 
procedure to obtain T'^{v) removes x and its adjacent edge in Tji{v). 



resp. black to those vertices of Tr{v) corresponding to elementary resp. non- 
elementary faces of R, see Figure El We identify each edge in Tr{v) with the 
corresponding edge in T(t>). 

To ease the presentation of our ideas, we assume for now that only cycles 
from TiiVj) are encountered in line 3 of the recursive greedy algorithm. In 
Section 15.41 we show how to handle cycles from B'^ U B'2 as well. 

So consider some iteration of the algorithm where a cycle C = C{v, e) G 
7i{Vj) has just been picked in line 3 and assume that all cycles added to B 
so far all belong to HiVj). Cycle C should be added to B only if ,B U {C} is 
nested. We will now show how to detect whether this is the case using the 
contracted dual trees. 

If there is a region R containing v such that Tr{v) contains e then e (in G) 
belongs to R (since otherwise, e would have been contracted in Tb,{v)). Since 
each cycle in B is isometric and since shortest paths are unique. Lemma [6] 
implies that S U {C} is nested. And the converse is also true: if i3 U {C} is 
nested then there is a region R containing C. In particular, R contains e so 
this edge must belong to Tji{v). 

It follows that detecting whether B U {C} is nested amounts to checking 
whether e is present in Tr{v) for some region R. 

Now, assume that B U {C} is nested (otherwise, we can discard C) and 
let us see how the contracted dual trees can help us check the condition in 



13 



line 4 of the recursive greedy algorithm. 

Define R to be the region containing C. Since e belongs to i?, this edge 
belongs to the contracted dual tree Tr{v). Let vi and V2 be the end vertices 
of e in Tpiiy). Removing e from Tji{v) splits this tree into two subtrees, one, 
Ti, attached to vi and one, T2, attached to V2- By Lemma [5l the condition 
in line 4 is satisfied if and only if Ti and T2 each contain at least one white 
vertex. 

Unfortunately, both of these two subtrees may contain many black ver- 
tices so for performance reasons, a simple search in these trees to determine 
whether they contain white vertices is infeasible. 

We therefore introduce pruned (contracted) dual tree T'j^{v), defined as 
the subtree of T^{v) obtained by removing a black degree one vertex and 
repeating this procedure on the resulting tree until all degree one vertices 
are white, see Figure [61 We refer to this as the pruning procedure. 

Lemma 7. With the above definitions, e G ) if and only if Ti and T2 
both contain white vertices. 

Proof. If Ti contains only black vertices then the pruning procedure will 
remove all vertices in Ti. In particular, the procedure removes Vi. Similarly, 
if T2 contains only black vertices then V2 is removed. In both cases, e is 
removed so e ^ T'^{v). 

Conversely, if both Ti and T2 contain white vertices then the pruning 
procedure does not remove all vertices from Ti and does not remove all 
vertices from T2. Hence, neither vi nor V2 is removed so e G T'j^{v). □ 

Lemma [7] shows that once T'^{v) is given, it is easy to determine whether 
both Ti and T2 contain white vertices and hence whether the condition in 
line 4 is satisfied: simply check whether e G T'pj{v). 

Note that if line 4 is satisfied, e G T'^{v) and hence e G Tr{v). By the 
above, this implies that B U {C} is nested. This shows that we only need 
T^{v) to test the condition in line 4. 

5.2 Inserting a Cycle 

In the previous section, we introduced contracted and pruned dual trees and 
showed how the latter can be used to test the condition in line 4 of the 
recursive greedy algorithm for cycles in Ti-lVj). In the following, we show 
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Figure 7: Faces of R belonging to Ri resp. R2 are identified by visiting the 
subtree of contracted dual tree Tr{v) consisting of gray resp. white vertices. 

how to maintain regions and contracted and pruned dual trees when such 
cycles are added to B in line 5. 

Initially, i3 = so the contracted and pruned dual trees are simply the 
dual trees T{v) for each boundary vertex v G Vj. And there is only one 
region, namely the external region R^{B). 

Now, suppose C = C{v,e) G H{Vj) has just been inserted into B in line 
5, see Figure [71 Let R be the region such that C splits R into internal region 
Ri and external region R2. We need to identify the faces of R belonging 
to Ri and to i?2- This can be done with two searches in contracted dual 
tree Tji{v). One search starts in the end vertex of e belonging to int{C) and 
avoids e (visiting the gray vertices in Figure [7j). The other search starts in 
the end vertex of e belonging to extiC) and also avoids e (visiting the white 
vertices in Figure [7]). It follows from Lemma [1] and from the definition of 
contracted dual trees that the first search identifies the faces of R that should 
belong to Ri and the second search identifies those that should belong to i?2- 

We also need to form one new face for namely the face defined by 
extiC). We denote this face by fu-^. Similarly, we need to form a new face 
for i?2, defined by int{C), and we denote this face by f^^. 

Next, we update contracted dual trees. The only ones affected are those of 
the form Tr{u), where u G R. There are three cases to consider: u G int{C), 
u G ext{C), and u E C. 

Case 1: Consider first a contracted dual tree Tf>{u) with u G intiC). Then 
u G -Ri so we need to discard Tf>{u) and construct Tn^iu). We obtain the 
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latter from the former by contracting all edges of Tr{u) having both end 
vertices in extiC) to a single vertex (this is possible by Lemma [1]). We 
identify this new vertex with the new face fn^ of Ri. 

Case 2: Now, assume that u G extiC). Then n G -R2 so Tji{u) should be 
replaced by Tr^{u). We do this by contracting all edges of Tr{u) having both 
end vertices in intiC) to a single vertex (again, we make use of Lemma [1]) 
and we identify this vertex with the new face f^^ of R2- 

Case 3: Finally, assume that u E C. Now, u belongs to both Ri and R2 
so we need to discard Tji{u) and construct Tr^{u) and Tji^{u). To do this, 
we first identify the edge e' in Tr{u) having one end vertex ui in int{C) and 
one end vertex U2 in ext{C). Then we construct the two trees Ti and T2 
formed by removing e' from Tr{u) with Ui G Ti and U2 G T2. We let T[ be 
Ti augmented with the edge from ui to fn-^ and let be T2 augmented with 
the edge from M2 to fn^ . 

It follows from Lemma [T] that T[ is the contracted dual tree Tr^ (m) for 
Ri and that is the contracted dual tree Tr^{u) for R2. 

We have described how to update contracted dual trees when C is added 
to B. We apply the same method to update pruned dual trees. The only 
difference is that the pruning procedure needs to be applied whenever a 
change is made to a pruned dual tree. 

5.3 Implementation 

Above, we gave an overall description of the algorithm when only cycles of 
TiiVj) are considered. We now go into more details and show how to give an 
efficient implementation of this algorithm. We start by describing the data 
structures that our algorithm makes use of. The main objects involved are 
regions, contracted dual trees, and pruned dual trees and we consider them 
in the following. 

5.3.1 Regions 

Associated with a region i? is a face list T{R) which is a linked list containing 
the faces of R. An entry of J^(-R) corresponding to a face / is assigned the 
colour white resp. black if / is elementary resp. non-elementary. If it is 
black, it has a bidirected pointer to the child of R contained in /. This gives 
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Data structure for region _R 




Data structure for Tf>{v) 



Figure 8: Illustration of data structures and some of their associated pointers. 

a representation of the region tree T{B). If the entry is white, it points to 
the corresponding elementary face of G. The entry also points to the entire 
data structure for R. 

Associated with the /-entry of J^{R) is also an array ^_r(/) with an entry 
for each boundary vertex in Vj. The entry of for a boundary vertex 

V belonging to R has a bidirected pointer to vertex / in contracted dual tree 
Tr(u), see Figure [HI It also has a bidirected pointer to vertex / in pruned 
dual tree T'^{v) if that vertex has not been deleted by the pruning procedure. 
All other entries of An^f) point to null. 

5.3.2 Contracted and pruned dual trees 

Associated with a contracted dual tree Tr{v) is a vertex list V{TR{y)) which 
is a linked list with an entry for each vertex of Tr{v). The entry for a vertex u 
points to the entry of T{R) for the face of R corresponding to u. Associated 
with the u-entry of V{Tji[v)) is also an edge adjacency list £u{Tii{v)), a linked 
list representing the edges adjacent to u in Tr{v). Each list entry contains a 
pointer to the u-entry of vertex list V(Tr(i;)) (allowing us to find the head 
of Su{T{i{v)) in constant time) as well as a bidirected pointer to an edge 
data structure. The edge data structure thus contains two pointers, one for 
each of its end vertices. Furthermore, it contains a bidirected pointer to the 
corresponding edge in dual tree T(w), see Figure [HI 

We keep a similar data structure for pruned dual tree T'i^{v). Both data 
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structures need to support edge contractions, edge insertions, and edge dele- 
tions and the data structure for T^{v) also needs to support the pruning 
procedure. We describe how to do this in the following. 

Edge contraction: We only describe edge contractions for contracted dual 
trees since pruned dual trees can be dealt with in a similar way. Assume we 
have a set Ec of edges (or edge data structures) in Tji{v) to be contracted 
to a single new vertex Vc and that these edges span a subtree of Tr(i;). We 
assume that we have a pointer to the entry of J^{R) corresponding to v^. 

To contract an edge e G E^, we first remove the pointer to the edge of 
dual tree T{v) corresponding to e. Traversing the two pointers associated 
with e, we find an entry in £ui(Tr{v)) and an entry in Su2(Tii{v)), where ui 
and U2 are the end vertices of e in TR(t>). 

We remove those two entries in lists Li = E^^ {T^{v)) and = £u2 {Tr{v)) 
and then merge Li and L2 to one list L since the new vertex is adjacent to 
edges adjacent to ui and U2 except e. If Li is appended to the tail of L2, 
we make every entry in Li point to the U2-entry in vertex list V(TR(f)). 
Otherwise, we make every entry in L2 point to the ui-entry in that list. For 
performance reasons, we append the shorter of the two lists to the tail of the 
other. 

We repeat the above for each edge of Ec and we end up with a single 
entry in V(Tf}(f)) representing the new vertex Vc- We make this entry point 
to the entry of J^{R) corresponding to Vc and we update the pointer to the 
f-entry in the associated array. 

How long does it take to contract edges? We will need the following 
lemma in our analysis (the proof can be found in the appendix). 

Lemma 8. Consider a set of objects, each assigned a positive integer weight. 
Let merge{o,o') be an operation that replaces two objects and 0' by a new 
object whose weight is the sum of the weights of o and d . Assume that the 
time to execute merge{o,o') is bounded by the smaller weight of objects 
and 0' . Then repeating the merge-operation on pairs of objects in any order 
until at most one object remains takes 0(14^ log W) time where W is the total 
weight of the original objects. 

Fix a. V & Vj and consider the set of contracted dual trees of the form 
Tr{v) generated during the course of the algorithm. Each time a cycle from 
T-iiVj) is added to B, at most two new edges are inserted into trees of this 
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form (case 3 in Section [32]) • Hence, there are 0{n) edges in total. It then fol- 
lows easily from Lemma [H] and from the way we concatenate lists during edge 
contractions that the total time spent on edge contractions in all contracted 
dual trees of the form Tii{v) is 0{nlogn). Since the number of choices of 
V is 0(\/n), we get a bound of 0(n'^/^ log n) time for all edge contractions 
performed by the algorithm. 

Edge deletion: We also describe this only for contracted dual trees. So 
suppose we are to delete an edge e = (ui,M2) from Tji{v). We need to form 
two new trees, Ti and T2. Let Ti be the tree containing ui and let T2 be the 
tree containing U2- For i = 1,2, a simple search (say, depth-first) in Tfi{v) 
starting in Ui and avoiding e finds the vertices of Tj in time proportional to the 
size of this tree. By alternating between these two searches (i.e., essentially 
performing them in parallel), we can find the vertices of the smaller of the 
two trees in time proportional to the size of that tree. 

Suppose that, say, Ti is the smaller tree. Then we can form the two data 
structures for Ti and T2 in time proportional to the size of Ti: extract the 
entries of vertex list V(Tr(w)) that should belong to Ti and form a new vertex 
list containing these entries. The old data structure for Tji{v) now becomes 
the new data structure for T2 after the entries have been removed. We also 
need to remove the pointer between e and the corresponding edge in dual 
tree T{v) and remove e from the edge adjacency lists but this can be done 
in constant time. 

The following lemma, which is similar to Lemma [HI immediately implies 
that the total time for edge deletions is O^n^^"^ logn) (the proof of the lemma 
is in the appendix). 

Lemma 9. Consider an object o with a positive integer weight W . Let split 
be an operation that splits an object of weight at least two into two new objects 
of positive integer weights such that the sum of weights of the two equals the 
weight of the original object. Assume that split runs in time proportional to 
the smaller weight of the two new objects. Then repeating the split -operation 
in any order, starting with object o, takes 0{W\ogW) time. 

Edge insertion: The only situation where edge insertions are needed is in 
case 3 of Section 15.21 With our data structure, this can clearly be done in 
constant time per insertion. 
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Pruning procedure: Finally, let us describe how to implement the prun- 
ing procedure for pruned dual trees. Recall that this procedure repeatedly 
removes black degree one vertices until no such vertices exist. 

We only need to apply the pruning procedure after an edge contraction 
and after an edge deletion (edge insertions are not needed in pruned dual 
trees since these edges will be removed by the pruning procedure). Let us 
only consider edge deletions since edge contractions are similar. 

Consider a pruned dual tree T'^{v) and suppose the algorithm removes 
an edge e = (mi,M2) from this tree. This forms two new trees Ti and T2, 
containing ui and M2, respectively. In Ti, only ui can be a black degree 
one vertex since in T^{v), no vertices had this property. Checking whether 
Ml should be removed takes constant time. If it is removed, we repeat the 
procedure on the vertex that was adjacent to ui. We apply the same strategy 
in T2, starting in U2. 

The total time spent in the pruning procedure is proportional to the 
number of vertices removed. Since the number of vertices only decreases and 
since the initial number of vertices in all pruned dual trees is 0{n^^'^), the 
total time spent by the pruning procedure is 0{n^^'^). 

5.3.3 The algorithm 

Having described the data structures involved and how they can support the 
basic operations that we need, let us show how to give an efficient imple- 
mentation of our algorithm. Still, we only consider cycles from T-C{Vj) in the 
for-loop. 

Initialization: First, we consider the initialization step. Applying the sep- 
arator theorem of Miller gives us and Vj in linear time. For each boundary 
vertex v, we need to compute shortest path tree T{v) and shortest path dis- 
tances from V in G. This can be done in 0{nlogn) time with Dijkstra's 
algorithm for a total of 0(?t,^/^ logn) time (in fact, a shortest path tree can 
be computed in linear time [13] but this will not improve the overall running 
time of our algorithm). We also need to compute dual trees T[v) and this 
can easily be done in 0(n^/^) additional time. These dual trees are also the 
initial contracted and pruned dual trees. Since we need all three types of 
trees during the course of the algorithm, three copies of each dual tree are 
initialized. 
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The algorithm then recursively computes Bi and 82- It is assumed that 
the recursive calls also return the weights of cycles in these sets. 

Our algorithm needs to extract B[ and B2 from these sets. This is done 
as follows. For every shortest path tree T{v) that has been computed in 
recursive calls (we assume that these trees are kept in memory), we mark 
vertices of T{v) belonging to Vj. Then we mark all descendants of these 
vertices in T{v) as well. Now, a Horton cycle C{v,e) obtained by adding e 
to T{v) contains a vertex of Vj if and only if at least one of the end vertices 
of e is marked. Since the total size of all recursively computed shortest path 
trees is bounded by the total space requirement which is 0(n^/^), it follows 
that B'l and B2 can be extracted from Bi and B2 in 0{n^^^) time. 

The cycles in U7Y(Vj) need to be sorted in order of non- decreasing 

weight. We are given the weights of cycles in B[ U B2 from the recursive calls 
and we can compute the weights of cycles in 7i(Vj) in a total of 0(n^/^) time 
using the shortest path distances computed above. Hence, sorting the cycles 
in B[U B'2U 'H{Vj) can be done in 0(?7,^/^ logn) time. 

Testing condition in line 4: Next, we consider the for-loop of the algo- 
rithm for some cycle C = C{v,e) E 'HiVj). As we saw in Section I^TTl testing 
the condition in line 4 amounts to testing whether dual edge e in T[v) is 
present in some pruned dual tree. Recall that we keep pointers between 
edges of dual trees and pruned dual trees. Since we remove a bidirected 
pointer between an edge data structure and the corresponding edge in a dual 
tree whenever it is contracted or deleted in a pruned dual tree, we can thus 
execute line 4 in constant time. 

Inserting a cycle: Line 5 requires more work and we deal with it in the 
following. Suppose we are about to add the above cycle C to i3 in line 5. 
With the pointer associated with e, we find the corresponding edge data 
structure in a contracted dual tree Tr(i;). Traversing pointers from this data 
structure, we find the data structure for R in constant time. This region 
should be split into two new regions Ri and R21 where Ri is the internal and 
R2 the external region w.r.t. R and C . We need to identify the boundary 
vertices and the of faces in R that belong to Ri and i?2, respectively. 

Identifying boundary vertices in Ri and R2: We first identify the set 
5{R) of boundary vertices of Vj belonging to R by traversing any one of 
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the arrays ArIJ) associated with an entry of T{R) and picking the vertices 
not having null-pointers. This takes 0{^/n) time. Since the total number of 
times we add a cycle to B is 0{n), total time for this during the course of 
the algorithm is 0(n'^/^). 

We will extract three subsets from 5{R): the subset 6int{R, C) of vertices 
belonging to mt{C), the subset 6ext{R, C) belonging to ext{C), and the subset 
5{R,C) belonging to C. 

If we can find these three subsets, we also obtain sets S{Ri) and 6{R2) 
of boundary vertices for i?i and R2, respectively, since S{Ri) = 6{R,C) U 
5,nt{R, C) and S{R2) = 5{R, C) U S,,t{R, C). 

The following lemma bounds the time to find the three subsets. The 
proof is somewhat long and can be found in the appendix. 

Lemma 10. With the above definitions, we can find in 0{^/n) time the sets 
5int{R,C), 6extiR,C), and 6{R,C) with 0{n^/'^ log n) time and 0{n^/'^) space 
for preprocessing. 

Lemma fTOl implies that the total time spent on computing sets of bound- 
ary vertices over all regions generated by the algorithm is Oin?/"^) (plus 
0{Tn?/'^\ogn) time for preprocessing). 

Identifying faces of Ri and R2: Having found the boundary vertices 
belonging to Ri and i?2, we next focus on the problem of identifying the 
faces of R belonging to each of the two new regions. 

As previously observed (see Figure [7]), we can identify the faces of Ri 
resp. R2 with, say, a depth-first search in Tjiiy) starting in the end vertex of 
e belonging to int{C) resp. ext{C) and avoiding e. We use the edge adjacency 
lists to do this. By alternating between the two searches, we can identify the 
smaller set of faces in time proportional to the size of this set. 

Let us assume that internal region Ri contains this smaller set (the case 
where external region R2 contains the set is similar). The search in Tfi{v) 
visited the entries of V(Tr(?;)) corresponding to faces in Ri. Since each such 
entry points to the corresponding entry in J^{R), we can thus identify the 
faces in this face list that should belong to J-'{Ri). 

We can extract these faces in time proportional to their number and thus 
form the face lists T{Ri) and J-'{R2) in this amount of time. By reusing the 
arrays associated with entries of J^{R), we do not need to form new arrays for 
jF(i?i) and jF(i?2)- However, we need to set the pointers of some entries of 
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these arrays to null. For the new null-pointers are those corresponding to 
boundary vertices of dext{R, C) since these are the boundary vertices of R not 
belonging to -Ri. And for R2, the new null-pointers are those corresponding 
to boundary vertices of Sint{R, C). 

Since we index the arrays by boundary vertices, we can identify pointers 
to be set to null in constant time per pointer. Pointers that are set to null 
remain in this state so we can charge this part of the algorithm's time to the 
total number of pointers which is 0{n^^'^). 

We also need to associate a new face with the data structure for Ri and 
for i?2 (i.e., faces //j^ and fn^ in Section [512]) . And we need to initialize an 
array for each of these two faces. This takes 0{y/n) time which is 0(n^/^) 
over all regions. 

Contracted and pruned dual trees for Ri and R2: What remains is to 
construct contracted and pruned dual trees for Ri and R2. Due to symmetry, 
we shall only consider contracted dual trees. We have already given an overall 
description of how to do this in Section 15. 2[ As we showed, 

1. for each u G Sint{R, C), we obtain Tr^^u) from Tr{u) by contracting all 
edges belonging to ext{C), 

2. for each u G 6ext{R, C), we obtain Tn^iu) from Tuiu) by contracting all 
edges belonging to int{C), and 

3. for each u G S{C), we obtain Tr^ (u) and (u) from Tr{u) by removing 
the unique edge in Tuiu) having one end vertex in int{C) and one end 
vertex in ext{C). 

In Section [5.3.2[ we described how to support edge contraction, edge deletion, 
and edge insertion such that the total time is 0(?7,^/^ logn). The only detail 
missing is how to efficiently find the edges to be contracted or removed in 
the three cases above. We consider these cases separately in the following. 

Case 1: Assume that Sint{R, C) ^ ^ and let u G Sint{R, C). 

With a depth-first search in Tr(w) as described above, we can identify all 
faces of R belonging to ext{C) in time proportional to the number of such 
faces. We can charge this time to the number of edges in Tr(m) that are to 
be contracted. 
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For each such face /, we can mark the corresponding vertex in T^{u) by 
traversing the pointer associated with entry u of array Anif)- Again, we can 
charge the time for this to the number of edges to be contracted. 

Now, we need to contract all edges of Tr{u) whose end vertices are both 
marked. In order to do this efficiently, we need to make a small modification 
to the contracted dual tree data structure in Section I5.3.2[ 

More precisely, we make the contracted dual trees rooted at some vertex. 
The choice of root is not important and may change during the course of 
the algorithm. What is important is that each non-root vertex now has a 
parent. By checking, for each marked non-root vertex whether its parent is 
also marked, we can identify the edges to be contracted in time proportional 
to the number of such edges. Of course this only works if the parent of a 
vertex can be obtained in constant time. Let us show how the contracted 
dual tree data structure can be adapted to support this. 

Recall that each vertex of a contracted dual tree Tr(v) is associated with 
an edge- adjacency list £u{Tji{v)) containing the edges adjacent to u in Tji{v). 
We now require the edge from v to its parent (if defined) to be the located 
at the first entry of this list. This allows us to find parents in constant time. 

How do we ensure that the parent edge is always located at the head of 
the list? This is not difficult after an edge insertion or deletion so let us focus 
on edge contractions. When an edge e = (mi,M2) is contracted, either ui is 
the parent of U2 or U2 is the parent of Ui. Assume, say, the former. Then the 
parent of Ui becomes the parent of the new vertex obtained by contracting 
e. When the two edge adjacency lists are merged, one of the two heads of 
the two old lists should thus be the head of the new list. This can easily be 
done in constant time. 

Case 2: This case is similar to case 1. 

Case 3: We need an efficient way of finding the unique edge e in Tr{u) 
having one end vertex in int{C) and one end vertex in ext{C). We do as 
follows: first we mark the entries in J^{R) corresponding to the set of faces 
of R belonging to int{C) or the set of faces of R belonging to ext{C). The 
set we choose to mark is the smaller of the two. We do this with "parallel" 
searches in Tr(w) as described above, using time proportional to the number 
of marked faces. 

We mark the corresponding vertices of Tr[u) (using pointers from the 
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arrays associated with entries of J-'{R)). By Lemma [H these form a subtree 
of Tji{u) so we can find e by starting a search in any marked vertex of Tji{u) 
and stopping once we encounter a vertex which is not marked. Then e is the 
last edge encountered in the search. This search also takes time proportional 
to the number of marked faces. 

Hence, constructing the contracted and pruned dual trees for Ri and 
R2 takes time proportional to the number of marked faces. Lemma M then 
implies that the total time for this during the course of the algorithm is 
0(n^/^ logn). 

Having constructed the contracted and pruned dual trees for Ri and R2, 
what remains before adding C to i3 is to add bidirected pointers between 
entries of the array associated with the new face in T{Ri) resp. T{^R2) and 
the new vertex in the contracted/pruned dual tree for Ri resp. R2- Since 
the size of the array is 0{\/n)^ this can clearly be done in a total of 0(?t,^/^) 
time. 

This concludes the description of the implementation of our algorithm. 
We have shown that it runs in 0(n^/^ logn) time and requires 0(n^/^) space. 

5.4 Recursively Computed Cycles 

So far, we have assumed that only cycles from 'HiVj) are encountered in line 
3 of the recursive greedy algorithm. Now, we show how to deal with cycles 
from B'l U B2. In the following, we only consider B[ since dealing with B2 is 
symmetric. 

The overall idea is the following. When a cycle C G Ti-iVj) is added to B, 
all cycles of B[ that cross C are marked. If in the for-loop, a cycle C & B[ is 
picked, it is skipped if it is marked since the GMCB is nested by Lemma O 
Otherwise, C must be fully contained in some region of the form R{C',B), 
C e B. Then C is added to B if and only if C separates a pair of elementary 
faces of R{C\B). 

We will assume that the recursive invocation of the algorithm in Gi re- 
turns region tree T{Bi) in addition to Bi. 

By applying Lemma [31 we see that every pair of elementary faces of Gi is 
separated by some cycle of Bi. Hence, each region associated with a vertex 
u of T{Bi) contains exactly one elementary face of Gi and we assume that 
the recursive call has associated this face with u. We let R{f, Bi) denote the 
region containing elementary face /. 
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(a) 



Figure 9: (a): Neither R{C',Bi), R{C",Bi), nor R{C"',Bi) are ancestors of 
R{fj, Bi) and only R{C', Bi) and R{C", Bi) are ancestors of both Bi) 
and Rif2,Bi). Thus, C crosses C" and C" and not C". (b): Both R{C',Bi) 
and R{C",Bi) are ancestors of R{fj,Bi) and only R{C',Bi) is an ancestor 
of neither -Bi) nor -R(/2, Si). Thus, C crosses C and not C". 

We use the conditions in the following lemma to identify those cycles of 
B[ that should be marked whenever a cycle of TiiVj) is added to B. 

Lemma 11. Let C = C{v,e) G 'HiVj). If e does not belong to Gi then C 
does not cross any cycle of B[. Otherwise, let fi and /2 be the elementary 
faces of Gi adjacent to e and let fj be the elementary face of Gi containing 
J . Then the set of cycles G' e B'^ that G crosses are 'precisely those which 
satisfy one of the following two conditions: 

1. R{G',Bi) is not an ancestor of R{fj,Bi) and is an ancestor of both 
R{fi,Bi) andR{f2,Bi) inT{Bi) (Figure\E(a)), 

2. R{G',B i) is an ancestor of R{fj,Bi) and is an ancestor of neither 
R{fi,B^) norR{f2,B,) znT{Bi) (Figure\^b)). 

The proof can be found in the appendix. 

The next lemma will simplify the test in line 4 of the recursive greedy 
algorithm for G & B'^. Again, the proof is in the appendix. 



26 



Lemma 12. Suppose that in the recursive greedy algorithm, C ^ B[ is the 
cycle currently considered and assume that it does not cross any cycle of the 
partially constructed GMCB B of G. If J C. ext{C) then all descendants of 
C in region tree T{Bi) belong to the GMCB of G. If J C int{C) then all 
cycles of non-descendants of G in T{Bi) belong to the GMGB ofG. 

Now, we are ready to describe how the algorithm deals with cycles from 
B'^. Each cycle in this set is in one of three states: active, passive, or cross 
state. 

Initially, all cycles in B[ are active. When a cycle from TiiVj) is added 
to B, Lemma [TT] is applied to identify all cycles from B'^ that cross this cycle. 
These cycles have their state set to the cross state. 

When the algorithm encounters a cycle C G in the for-loop, G is 
skipped if it is in the cross state. 

If G is active, it is completely contained in some region R. There are two 
cases to consider: J C ext{C) and J C int{C). We assume that J C ext{G) 
since the case J C int{C) is similar. We need to determine whether G should 
be added to B. By Lemma [5|, this amounts to checking whether there are two 
elementary faces of R which are separated by G . By Lemma [121 we know 
that the elementary faces of R belonging to int{G) are exactly the elementary 
faces of the region R' in intiC) that was generated when G was added to Bi 
during the recursive call for Gi. 

Hence, we add C to S if and only if the number of elementary faces in R 
is strictly larger than the number of elementary faces in R' . 

If G is added to i3, region R is split into two smaller regions. Let -Ri be 
the internal region and let R2 be the external region. Since J C ext{G), 
Lemma [12] implies that the cycles belonging to int{C) that are added to 
B during the course of the algorithm are exactly G and its descendants in 
T{Bi). We therefore do not need to maintain Ri or any regions contained in 
mt{G). 

Instead, we make all descendants of G in T{Bi) passive. When a passive 
cycle is encountered by the algorithm, there is no need to update regions or 
contracted or pruned dual trees and the cycle is simply added to B. 

Now, let us consider R2. In order to obtain this region, we replace all 
faces of R belonging to int{C) with a single new face defined by intiC). And 
we contract all edges in int{G) to a single black vertex in all contracted and 
pruned dual trees for R. 

This completes the description of the extension of our algorithm that 
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deals with cycles from B[U 82- 
5.4.1 Implementation 

Let us show how to give an efficient implementation of the above algorithm 
for cycles from B'lU 82- Due to symmetry, we may restrict our attention to 
B[ in the following. 

Identifying cross state cycles: The ffist problem is to identify the cycles 
of B[ that should be in the cross state when a cycle C = C{v, e) G TiiVj) is 
added to B. 

To solve this problem, we apply Lemma [TTJ Checking whether e belongs 
to Gi takes constant time. If e is not an edge of Gi then no new cycles will 
be in the cross state. Otherwise, we obtain elementary faces /i and /2 in 
constant time since these are the end vertices of e in the dual of Gi. 

We assume that we can compute lowest common ancestors in T{Bi) effi- 
ciently. We can use the data structure of Harel and Tarjan [TT] for this. 

Let ai be the lowest common ancestor of R{fi,Bi) and R{f2,Bi) in 
T{Bi), see Figure [TOl Let 02 be the lowest common ancestor of R{fi, Bi) and 
R{fj, Bi). Let 03 be the lowest common ancestor of R{f2, Bi) and R{fj, Bi). 
Finally, let P be the path in T{Bi) containing R{fj,Bi) and its ancestors. 

A cycle C G B[ satisfies the ffist condition in Lemma [TT] if and only if 
it is not associated with a vertex on P and if it is associated with ai or an 
ancestor of ai (Figure [TOT a)). And it satisfies the second condition if and 
only if it is associated with a vertex on P and not with a2, 03, or an ancestor 
of either of these two vertices (Figure [TOlfb)). 

To identify cycles that satisfy the ffist condition, we start at ai and walk 
upwards in T{Bi), marking cycles as we go along. The process stops when a 
vertex on P is reached. 

To identify cycles satisfying the second condition, we instead move up- 
wards in T{Bi) along P, starting in R{fj,Bi). We stop when the root of 
T{Bi) or when 02 or 03 is reached. 

Although this strategy works, it is slow since the same cycles may be con- 
sidered several times during the algorithm. To remedy this, we ffist observe 
that when identifying cycles associated with vertices from ai to P, we may 
stop if we encounter a cycle that is already in the cross state since then all 
its ancestors which are not on P must also be in this state. 
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Figure 10: (a): Cycles associated with ai or an ancestor of ai and not with 
a vertex on P are exactly those that satisfy the first condition in Lemma [TTJ 
(b): Cycles associated with a vertex on P and not with 02, 03, or an an- 
cestor of either a2 or are exactly those satisfying the second condition in 
Lemma [TTJ 

Next, we observe that when identifying cycles associated with vertices on 
P, we always consider them from bottom to top. Hence, by keeping track 
of the bottommost b vertex on P whose associated cycle is not in the cross 
state, we can start the next traversal of P from b. If the cycle associated with 
02 or with 03 is already in the cross state, we need not consider any vertices. 
Otherwise, we walk upwards in P from b, changing the state of cycles to the 
cross state and stop if 02 or 03 is reached. 

It follows that we can identify cycles satisfying one of the two conditions 
and change their state in time proportional to the number of cycles whose 
state chang result of this. Hence, the total time for this is bounded 

by the size of T{Bi) which is linear. 

Testing condition in line 4: In the following, let C be an active or passive 
cycle in B[ just encountered by our algorithm. We will assume that J' C 
ext{C). The case JT" C intiC) is similar. 

We first need to determine whether C should be added to B. This is trivial 
if C is passive since passive cycles should always be added. And as noted in 
Section EH no pruned dual trees need to be updated after the insertion of a 
passive cycle. 

So assume that C is active. Let R be the region containing C and let R' 
be the region in int{C) that was generated when C was added to Bi during 
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the construction of the GMCB of Gi. As we showed above, determining 
whether C should be added to B amounts to checking whether the number 
of elementary faces in R is strictly larger than the number of elementary 
faces in R' . 

We can easily extend our region data structure to keep track of the num- 
ber of elementary faces in each region without increasing the time and space 
bounds of our algorithm. By recording this information for R! during the re- 
cursive call for Gi, it follows that we can determine in constant time whether 
R contains more elementary faces than R' . 

Of course, this only works if we can quickly identify R and R' . Identifying 
R' is simple since this region is associated with the vertex vc of region tree 
T{Bi) associated with C. 

To identify i?, let R^^ be the region associated with vc in T{Bi). Since 
Bi is the GMCB of Gi, each pair of elementary faces of Gi is separated by 
some cycle of Bi. It follows that R^^ contains exactly one elementary face 

of Gi- We may assume that this face was associated with vc during the 
construction of Bi so that we can obtain this face in constant time from vc- 

Face is also an elementary face in G and it belongs to R. Recall 
from Section 15.3.11 that there is a bidirected pointer between R and each 
elementary face of G belonging to R. Hence, we can obtain R from /^^ in 
constant time 

It follows from the above that we can check if C should be added to B in 
constant time. 

Inserting a cycle: Now, suppose C should be inserted into B. We first 
make cycles of B'l passive according to Lemma [121 This can be done with, 
say, a depth-first search in T{Bi) starting in vertex vc of T{Bi) and visiting 
descendants of this vertex. The search stops when a vertex associated with a 
passive cycle is encountered. Each search identifies the vertices of T{Bi) that 
are associated with cycles whose state changes from non-passive to passive. 
And since we stop a search when a passive cycle is encountered, all searches 
take total time proportional to the size of T{Bi) which is 0{n). 

Next, we need to update regions and contracted and pruned dual trees. 
Let Ri be the internal region and let R2 be the external region w.r.t. R 
and C. As we showed in the overall description of the algorithm, the only 
problem that we need to consider is how to construct R2 and its contracted 
and pruned dual trees. We showed that this amounts to replacing all faces 
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Figure 11: Solid cycles belong to the partially constructed GMCB B of G. 
All descendants of vertices of T{Bi) associated with solid cycles are passive 
and thus do not define faces of R. For this instance, Ci and C2 define the 
non-elementary faces of R. 

of R belonging to int{C) with a single new face defined by int{C) and to 
contract all edges in int{C) to a single black vertex in all contracted and 
pruned dual trees for R. 

We will show how to find the faces of R belonging to intiC) in time 
proportional to their number. Applying the charging schemes introduced in 
Section I5.3.3[ this will suffice to prove the desired time and space bounds for 
the entire algorithm. 

Consider an active cycle C associated with a descendant u of vc in T{Bi). 
If C was previously considered in the for-loop of our algorithm, it must have 
been added to B (by Lemma [T2|) . This implies that int{C') must be a non- 
elementary face of R since otherwise, C would be passive, see Figure [TTl The 
converse holds as well: any non-elementary face of R belonging to intiC) is 
realized by int{C') for such a cycle C previously considered in the for-loop. 

It follows that we can find all non-elementary faces of R belonging to 
int{C) by identifying the active descendants of C in T{Bi) that have already 
been considered in the for-loop. Since all active cycles associated with de- 
scendants of C are to become passive, we can charge the time for finding these 
faces to the number of cycles whose state changes from active to passive. 

What remains is to find the elementary faces of R belonging to int{C). 
Recall that we have associated with each vertex u of T{Bi) the unique el- 
ementary face of Gi contained in the region associated with u. Vertex vc 
and its descendants in T{Bi) are thus associated with exactly the elementary 
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faces of Gi belonging to int{C). These faces are also elementary faces in G. 

It follows that the elementary faces of R belonging to int{G) are asso- 
ciated with exactly the descendants of G corresponding to active cycles not 
already considered by the algorithm. Using the same charging scheme as 
above, we can also identify these faces within the required time and space 
bounds. 

We have shown how to implement the entire recursive greedy algorithm to 
run in 0(n^/^ log n) time and 0(n^'^^) space and we can conclude this section 
with the main result of our paper. 

Theorem 2. Given an n-vertex planar, undirected graph G = {V,E) with 
non-negative edge weights, the following implicit representation of the GMCB 
B of G can be computed in 0{n^^'^logn) time and 0(n^^'^) space: 

1. a set of trees Ti, . . . , in G rooted at vertices Vi, . . . ,Vk, respectively, 

2. a set of triples {i, e, w) representing the cycles in B, where «€{!,..., k}, 
e = {u,v) E E \ Exi, and w G M, where u and v are vertices in Ti. The 
pair {i, e) represents the cycle in B formed by concatenating e and the 
two paths in Ti from Vi to u and from Vi to v, respectively. The value 
of w is the weight of this cycle, 

3. the region tree T{B) where each vertex points to the associated region, 

4- a set of regions. Each region is associated with the unique elementary 
face of G contained in that region. Each internal region R{G, B) is 
associated with the triple representing G . 

6 Corollaries 

In this section, we present results all of which follow from Theorem [2j The 
first is an output-sensitive sensitive algorithm for computing an MCB. 

Corollary 1. A minimum cycle basis of an n-vertex planar, undirected graph 
with non-negative edge weights can be computed in 0{n^^'^\ogn + G) time and 
Q 1^^3/2 _j_ space, where G is the total length of all cycles in the basis. 

Proof. Follows immediately from Theorem [2l □ 

A stronger result holds when the graph is unweighted: 
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Corollary 2. A minimum cycle basis of an n-vertex planar undirected, un- 
weighted graph can be computed in 0(n^/^Iogn) time and 0{n^^'^) space. 

Proof. Let G be an n-vertex planar, undirected, unweighted graph. The 
internal elementary faces of G define a cycle basis of of G of total length 
0{n). Hence, since G is unweighted, an MCB of G has total length 0{n). 
The result now follows from Corollary [1] □ 

Since the all-pairs min cut problem is dual equivalent to the MCB problem 
for planar graphs, we also get the following two results. 

Corollary 3. All-pairs min cuts of an n-vertex planar, undirected graph 
with non-negative edge weights can be computed in 0(n^/^logn + G) time 
and 0{n^''^ + G) space, where G is the total length of the cuts. 

Proof. Let G be an n-vertex planar, undirected graph with non-negative edge 
weights. As shown in [12], if G is connected, we can solve the APMCP for 
G by solving the MCBP for the dual G* of G. 

We may assume that G is connected since otherwise, we can consider each 
connected component separately. We cannot immediately solve the MCBP 
for G* since this is a multigraph. But we can avoid an edge of the form {u, u) 
by splitting it into two edges (m, f) and {v,u) whose sum of weights equal 
the weight of {u,u). And we can avoid multi-edges in a similar way. Let 
G' be the resulting planar graph. It is easy to see that G' has size 0{n). 
Furthermore, an MCB B of G' can be transformed into an MCB of G* in 
time proportional to the total size of cycles in B. The result now follows 
from Corollary [H □ 

Corollary 4. All-pairs min cuts of an n-vertex planar, undirected, unweighted 
graph can be computed in 0(n^/^logn) time and 0{n^^'^) space. 

Proof. This result is easily obtained by combining the proofs of Corollary [2] 
and Corollary [31 □ 

Next, we present our subquadratic time and space algorithm for finding 
the weight vector of a planar graph. 

Corollary 5. The weight vector of an n-vertex planar, undirected graph with 
non-negative edge weights can be computed in 0(n^/^ log n) time and Oin?/"^) 
space. 
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Proof. From Theorem [21 we obtain an implicit representation of the GMCB 
B for the input graph. We then compute the weights of all cycles in B using 
linear time and space. Sorting them takes 0(nlog?T.) time. This gives the 
weight vector of the input graph in a total of 0{n^^^\ogn) time and 0{n^^'^) 
space. □ 

From Theorem [2|, we also obtain a faster algorithm for computing a 
Gomory-Hu tree of a planar graph. 

Corollary 6. A Gomory-Hu tree of ann-vertex connected, planar, undirected 
graph with non-negative edge weights can be computed in 0{n^^'^ logn) time 
and 0{n^^'^) space. 

Proof. The following algorithm constructs a Gomory-Hu tree for G [26]: a 
tree T spanning a collection of vertex sets Si, . . . , St is maintained, starting 
with Si = V. At each step, a set Si is picked such that 15^1 > 1 and any 
two distinct vertices u,v E Si are chosen. Set Si is then regarded as the root 
of T and each subtree of T, i.e., each tree in T \ {Si}, is collapsed into a 
single supernode. A min u-v cut in the resulting graph is found, partitioning 
V into two subsets, Vi and V2, where u E Vi and v G V2. Tree T is then 
modified by splitting Si into two vertices, S'j^ and S'jj , where 5"^^ = Sir\Vi and 
Si^ = Si nV2. The two vertices are connected by a new edge whose weight 
equals the size of the min cut found. Finally, each subtree of the old T is 
connected to 5",^ if the corresponding supernode was in the same partition as 
u in the cut. Otherwise, the subtree is connected to Si^. 

Let us show how to implement this algorithm to obtain the desired time 
and space bounds. We first apply Theorem [2] to the dual G* of G, giving 
an implicit representation of the GMCB of G*. By Lemma [31 each cycle G 
in this basis is a minimum-weight cycle that separates some pair of faces fi 
and /2 in G*. Let Ui and U2 be the vertices of G corresponding to /i and /2, 
respectively. By duality of the GMCBP and the APMCP [12], the edges of 
G are the edges of a min Ui-Vi cut in G of weight equal to the weight of G. 

Now, pick any G cycle in the GMCB of G*. As the initial min cut in 
the Gomory-Hu tree algorithm, we pick the one corresponding to G. This 
separates the initial set Si = Si = V into two sets Si^ and Si^, where Si-^ is 
the set of vertices of G corresponding to faces of G* in int{C) and S'jj is the 
set of vertices of G corresponding to faces of G* in ext{G). Now, T consists 
of vertices S'j^ and Si^ and an edge (Sj^, SjJ. The weight of this edge is equal 
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to the weight of C. Since we are given the weight of C from Theorem [21 we 
can this obtain the weight of edge (S'j^, Si^) in constant time. 

Note that for each pair of vertices u and v in S'jj, there is a min u-v cut 
defined by a cycle of B which is a descendant of C in T{E). And for each 
pair of vertices u and v in Si^ , there is a min u-v cut defined by a cycle of ;B 
which is a non- descendant of C in T{B). 

Hence, we have separated our problem in two and we can recursively 
compute the Gomory-Hu tree for G by splitting region tree T{B) in two at 
each recursive step. The recursion stops when we obtain a set Si of size one. 
At this point, we obtain the elementary face / of G* correponding to the 
vertex in Si using part four of Theorem [2J The vertex of G corresponding to 
/ in G* is then the unique vertex in Si. 

Let us analyze the running time of this algorithm. Applying Theorem [2] 
takes 0(?7,^/^ logn) time and 0(n'^/^) space. Note that in the algorithm above, 
we do not need to compute the vertices in the 5'j-sets until they have size one. 
So each step of the algorithm, where the current S'j-set has size greater than 
one, can be implemented to run in constant time. And we can also execute 
a step where l^jl = 1 in constant time using the fourth part of Theorem [2] to 
find the vertex in Si. 

Since the GMCB of G* contains 0(n) cycles, it follows that the algorithm 
runs in linear time and space, in addition to the time and space in Theorem [2j 

□ 

Finally, we present our oracle for answering min cut queries. 

Corollary 7. LetG he an n-vertex planar, undirected graph with non-negative 
edge weights. With 0(n^/^logn) time and 0(?t,^/^) space for preprocessing, 
the weight of a min cut between any two given vertices of G can he reported 
in constant time. The cut itself can he reported in time proportional to its 
size. 

Proof. We may assume that G is connected since otherwise, we can consider 
each connected component separately. We first construct a Gomory-Hu tree 
T of G. By Corollary [HI this can be done in 0{rt'^'^\ogn) time and 0(n^/^) 
space. By definition of Gomory-Hu trees, answering the query for the weight 
of a min cut between two vertices u and v oi G reduces to answering the 
query for the minimum weight of an edge on the simple path between u and 
v in T. 
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It is well-known that any tree with m vertices has a vertex v such that 
the tree can be split into two subtrees, each rooted at v and each containing 
between m/4 and 3m/4 vertices. Furthermore, this separator can be found 
in linear time. 

We find such a separator in T and recurse on the two subtrees. We stop 
the recursion at level \og{^/n). The total time for this is 0{n\ogn). 

Let S be the subtrees at level \og{y/n). We observe that these trees are 
edge-disjoint and their union is T. Furthermore, |iS| = 0{^/n) and each 
subtree has size 0{y/n). The boundary vertices of a subtree S E S are the 
vertices that S shares with other subtrees in S. Vertices of S that are not 
boundary vertices are called interior vertices of S. We let B be the set of 
boundary vertices over all subtrees in S. It is easy to see that \B\ = 0{^/n). 

For each boundary vertex 6 e S, we associate an array with an entry for 
each vertex of T. The entry corresponding to a vertex v ^ b contains the 
edge of minimum weight on the simple path in T between b and v. 

Since = 0{y/n), it follows easily that we can construct all these arrays 
and fill in their entries in a total of 0(n^/^) time and space. This allows us 
to answer queries in T in constant time when one of the two vertices belongs 
to B. 

We associate each vertex w of T not belonging to B with the unique 
subtree in S containing v as an interior vertex. 

Associated with v is also an array with an entry for each S E S \ {S^}- 
This entry contains the vertex b of B belonging to such that any path 
from f to S* contains b. Note that for any other vertex v' of S^, any path 
from v' to S also contains b. From this observation and from the fact that 
\S\ = 0{y/n), it follows that we can compute the arrays associated with 
interior vertices in all subtrees using a total of 0{n^^^) time and space. 

Finally, we associate with v an array with an entry for each vertex v' of 
Si,. This entry contains the edge of minimum weight on the simple path in 
Sy from V to v'. Since has size 0{^/n), the entries in this array can be 
computed in 0{^/n) time. Over all interior vertices of all subtrees of S, this 
is 0(n^/^) time. 

Now, let us describe how to answer a query for vertices u and v in T. In 
constant time, we find the subtrees Su, Sy E S such that u E Su and v G Sy. 
If Su = Sy or if u or belongs to i?, we can answer the query in constant 
time with the above precomputations. 

Now, assume that 5*^ ^ Sy and that u and v are interior vertices. We 
find the boundary vertex b of S^ such that any path from u to Ry contains 
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b. Let Pi be the simple path in from m to & and let P2 be the simple path 
in T from b to v. For i = 1,2, the above precomputations allow us to find 
the least-weight edge Cj on Pi in constant time. Let e be the edge of smaller 
weight among ei and 62- Returning the weight of e then answers the query 
in constant time. 

To show the last part of the corollary, observe that when the weight of 
edge e is output by the above algorithm, the set of edges in the corresponding 
cut is defined by a cycle Ce in the GMCB B oi G*. During the construction 
of Gomory-Hu tree T (see Corollary [6]) , we can associate e with the implicit 
representation of Ce from Theorem [2l Hence, given e, we can output Cg in 
time proportional to its size. This completes the proof. □ 

7 Obtaining Lex-Shortest Path Trees 

Let w : E ^ M. he the weight function on the edges of G. In Section HI we 
assumed uniqueness of shortest path in G between any two vertices w.r.t. 
w. We now show how to avoid this assumption. We assume in the following 
that the vertices of G are given indices from 1 to n. 

By results in [12] , there is another weight function w' on the edges of G 
such that for any pair of vertices in G, there is a unique shortest path between 
them w.r.t. w' and this path is also a shortest path w.r.t. w. Furthermore, for 
two paths P and P' between the same pair of vertices in G, w'{P) < w'{P) 
exactly when one of the following three conditions is satisfied: 

1. P is strictly shorter than P' w.r.t. w, 

2. P and P' have the same weight w.r.t. w and P contains fewer edges 
than P', 

3. P and P' have the same weight w.r.t. w and the same number of edges 
and the smallest index of vertices in Vp\Vp' is smaller than the smallest 
index of vertices in Vp/ \ Vp. 

A shortest path w.r.t. w' is called a lex-shortest path and a shortest path tree 
w.r.t. w' is called a lex-shortest path tree. 

As shown in [12], lex-shortest paths between all pairs of vertices in G can 
be computed in 0{n^ logn) time. We need something faster. In the following, 
we show a stronger result, namely how to compute a lex-shortest path tree in 
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0(n log n) time. Since we only need to compute shortest paths from 0{y/ri) 
boundary vertices, this will give a total time bound of 0(n^/^ logn). 

We also need to find lex-shortest path trees in subgraphs of G when 
recursing and we need to compute them w.r.t. the same weight function w' . 
By the above, this can be achieved simply by keeping the same indices for 
vertices in all recursive calls. 

Now, let s G y be given and let us show how to compute the lex-shortest 
path tree in G with source s in O(nlogn) time. 

We first use a small trick from [12]: for function w, a sufficiently small 
e > is added to the weight of every edge. This allows us to disregard the 
second condition above. When comparing weights of paths, we may treat e 
symbolically so we do not need to worry about precision issues. 

We will apply Dijkstra's algorithm with a few additions which we describe 
in the following. We keep a queue of distance estimates w.r.t. w as in the 
standard implementation. Now, consider any point in the algorithm. Let 
d be the distance estimate function. Consider an unvisited vertex v with 
current distance estimate d[v\ < oo and predecessor vertex p. 

Suppose that at this point, the algorithm extracts a vertex p' from Q 
which is adjacent to w in G and suppose that d[p'] + w{p',v) = d[v]. The 
central problem is to decide whether v should keep p or get p' as its new 
predecessor. In the following, we show how to decide this in O(logn) time. 
This will suffice to give an 0(n log n) time algorithm that computes the lex- 
shortest path tree in G with source s. 

Let P be the path in the partially constructed lex-shortest path tree T 
from s to p followed by edge {p,v). Let P' be the path in T from s to p' 
followed by edge {p', v). Note that P and P' both have weight d[v] w.r.t. w. 
Hence, P is shorter than P' w.r.t. w' if and only if the third condition above 
is satisfied. In other words, v should keep p as its predecessor if and only if 
this condition is satisfied. 

Let q be the lowest common ancestor of p and p' in T. Paths P and P' 
share vertices from s to q. Then they split up and do not meet before v. 

Let Q be the subpath of P from the successor of q to p. Let Q' be the 
subpath of P' from the successor of q to p'. Testing the third condition above 
is equivalent to deciding whether the smallest vertex index in Vq is smaller 
than the smallest vertex index in Vq>. 

We assume that for each vertex u in T, we have pointers po[u], . . . ,Pky,[u] 
and values mo[u], . . . , mk^lu] G {!,..., n}. For i = 0, . . . , ku, Pi[u] points to 
the ancestor a of u in T for which the number of edges from a to u is 2*. 
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And mi[u] is the smallest vertex index on the path in T from a to u. The 
value of ku is defined as the largest i such that Pi[u] is defined. Note that 
ku = O(logn). 

Since P and P' have the same number of edges, the same holds for Q 
and Q'. From this observation, it follows that we can apply binary search 
on the pointers defined above to find lowest common ancestor q in 0(log?T,) 
time. And with these pointers and the mj-values, we can partition Q and Q' 
into O(logn) intervals in O(logn) time and find the smallest index in each 
interval in constant time per interval. Hence, we can decide whether the 
smallest vertex index in Vq is smaller than the smallest vertex index in Vqi 
in logarithmic time, which gives the desired. 

The only problem that remains is how to compute pointers and mj-values 
during the course of the algorithm. Whenever the partially constructed lex- 
shortest path tree is extended with a new vertex u, we need to compute 
Po[u], . . . ,Pku[u] and mofw], . . . ,mk^[u] G {!,... ,n}. But this can easily be 
done in O(logn) time using the pj-pointers and mj-values for the ancestors 
of u in T. 

We can now conclude this section with the following theorem. Since we 
did not make use of planarity in this section, we get a more general result, 
which we believe to be of independent interest. 

Theorem 3. A lex-shortest path tree in an undirected graph with m edges 
and n vertices can he computed in 0{{m + n) logn) time. 

Proof. Follows by combining the above with a standard implementation of 
Dijkstra's algorithm. □ 

8 Concluding Remarks 

We showed that finding a minimum cycle basis of an n-vertex planar, undi- 
rected, connected graph with non-negative edge weights requires fi(n^) time, 
implying that a recent algorithm by Amaldi et al. is optimal. We then pre- 
sented an algorithm with 0(n'^/^ log n) time and 0{rc'/'^) space requirement 
that computes such a basis implicitly. 

From this result, we obtained an output-sensitive algorithm requiring 
0{n^^'^\ogn + C) time and 0{rt'^'^ + C) space, where C is the total length 
of cycles in the basis that the algorithm outputs. For unweighted graphs, we 
obtained 0{rc^^'^\ogn) time and 0{n^/'^) space bounds. 
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Similar results were obtained for the all-pairs min cut problem for planar 
graphs since for planar graphs, this problem is known to be dual equivalent 
to the minimum cycle basis problem. 

As corollaries, we obtained algorithms that compute the weight vector 
and a Gomory-Hu tree of a planar n- vertex graph in 0(n^/^logn) time and 
0(?7,^/^) space. The previous best bound was quadratic. 

From the Gomory-Hu tree algorithm, we derived an oracle for answering 
queries for the weight of a min cut between any two given vertices. Pre- 
processing time is 0(n'^/^logn) and space is 0(n^/^). Quadratic time and 
space was previously the best bound for constructing such an oracle. Our 
algorithm can output the actual cut in time proportional to its size. 
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Appendix 



Proof of Lemmas [8] and [9] 

Let us first prove Lemma [H We only need to consider the hard case where in 
beginning, all objects have weight 1 and at termination, exactly one object 
of weight W remains. 

Consider running the algorithm backwards: starting with one object of 
weight W, repeatedly apply an operation split that splits an object of weight 
at least two into two new objects of positive integer weights such that the 
sum of weights of the two equals the weight of the original object. Assume 
that split runs in time proportional to the smaller weight of the two new 
objects. If we can give a bound of 0{W\ogW) for this algorithm, we also 
get a bound on the algorithm stated in the theorem. 

The running time for the new algorithm satisfies: 

T(w) < max iT(w') + T(w — w') + cw'} 

l<-i«'<L«>/2j 

for integer w > 1 and constant c > 0. It is easy to see that the right-hand 
side is maximized when w' = \_w/2\. This gives T(W) = 0(WlogW), as 
desired. 

The above proof also holds for Lemma M 
Proof of Lemma [TO] 

We need to show that for a cycle C = C{v,e) E TiiVj) belonging to a region 
-R, sets 6int{R,C), Sext{R,C), and S{R,C) can be computed in 0{y/n) time 
with 0(n^^^logn) preprocessing time and 0(n^/^) space. 

First, observe that since C is completely contained in R, 6{R, C) is the 
subset of all boundary vertices belonging to C. Hence, this subset does not 
depend on R. We will thus refer to it as 5{C) in the following. 

Let Wo, ... , Vr-i be the boundary vertices encountered in that order in a 
simple, say clockwise, walk of J and let J = Ji^Jx ■ ■ ■ Jr-i be a decompo- 
sition of J into smaller curves where Ji starts in Vi and ends in ^^(i+i) modr, 
z = 0, . . . , r — 1. Each curve J7i is completely contained in an elementary face 
of G and we let f{Ji) denote this face. 

In our proof, we need the following lemma and its corollary. 
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Lemma 13. Let P be a shortest path in G from a vertex u to a vertex v. 
Then a vertex w belongs to P if and only if dciu, w) + dc{w, v) = dc^u, v) . 

Proof. If w belongs to P then clearly dciu.w) + dciw.v) = dciu^v). And 
the converse is also true since shortest paths in G are unique. □ 

Corollary 8. Let C = C(f,e) be defined as above. Let w & V and assume 
that single-source shortest path distances in G with sources v and w have 
been precomputed. Then determining whether w belongs to C can be done in 
constant time. 

Proof. Let ui and U2 be the end vertices of e and let Pi resp. P2 be the 
shortest paths in G from v to ui resp. U2- Since G is isometric, both Pi 
and P2 belong to G and the union of their vertices is exactly the vertices of 
G. Hence, determining whether w belongs to G is equivalent to determining 
whether w belongs to Pi or to P2. The result now follows from Lemma [131 D 

We will assume that single-source shortest path distances in G with each 
boundary vertex as source have been precomputed. As observed earlier, this 
can be done in 0(ra^/^ logn) time and 0(n^/^) space. Corollary [8] then allows 
us to find the set S{G) of boundary vertices belonging to G in 0(r) = 0{y/n) 
time. We may assume that we have the boundary vertices on G cyclically 
ordered according to how they occur on in a clockwise walk of that curve. 

In the following, let Vi = v {so G = G{vi,e)). Consider two consecutive 
vertices Vi-^ and Vi^ of S{G) in this cyclic ordering. We assume that 12 i 
since the case 12= ^ can be handled in a similar way. There are two possible 
cases: 

1. the boundary vertices (excluding and Wjj) encountered when walking 
from Vi-^ to along J all belong to int{G), or 

2. they all belong to ext{C). 

Let fjg be the predecessor boundary vertex of Vi^ on J (i.e., is = {12 — 
1) mod r), see Figure [T2l Then elementary face f[Ji^) belongs to int{C) if 
and only if the first case above holds. This follows from the fact that J does 
not cross any edges of G. 

Lemma [T3] below shows how we can check whether fiJi^^) belongs to 
int{G). First, let u and v be the end vertices of e and let P„ and P^ be 
the shortest paths from Vi to u and v, respectively. Suppose w.l.o.g. that Vi^ 
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Figure 12: (a) The first condition and (b): the second condition in Lemma [T^ 



belongs to Pu, see Figure [121 Let u' be the predecessor of Vi^ on This 
is well-defined since 22 7^ ^- If 7^ u, let u" be the successor of Vi^ on P„. 
Otherwise, let u" = v (so m" is the vertex 7^ m' adjacent to on C). Let 
v' resp. f" be the predecessor resp. successor of in a clockwise walk of 

For three points p,q,r in the plane, let W{p,q,r) be the wedge-shaped 
region with legs emanating from p and with right resp. left leg containing q 
resp. r. 

Lemma 14. With the above definitions, f{Ji^) belongs to int{C) if and only 
if one of the following conditions hold: 

1. Pu is part of a clockwise walk of C (when directed from Vi to u) and 
W{vi^,u' ,u") contains W{vi^,v',v") (FigurellWa)). 

2. Pu is part of a counter-clockwise walk of C (when directed from Vi to 
u) and W{vi2,u" ,u') contains W{vi^,v' ,v") (FigureUWb)). 

Proof. Assume first that P„ is part of a clockwise walk of C, see Figure [T2](a). 
Then int{C) is to the right of the directed path u' Vi^ u". Since G is 
straight-line embedded, f{Jh^) belongs to int{C) if and only if W{vi^, u', u") 
contains W{yi.^,v' ,v"). 

Now, assume that Pu is part of a counter-clockwise walk of C, see Fig- 
ure [T2r b). Then int{C) is to the right of the directed path u" — > 
u'. Thus, f{Ji.^) belongs to int{C) if and only if W^Vi^^u" ,u') contains 
W{vi^,v',v"). □ 
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Lemma [T3] and the above discussion show that to efficiently determine 
whether the boundary vertices between and belong to int{C) or to 
ext{C), we need to quickly find u', u", v', and v" and determine whether 
is part of a clockwise or counter-clockwise walk of C. 

By keeping a clockwise ordering of vertices of all elementary faces, we can 
find v' and v" in constant time. For each shortest path tree in G rooted at a 
boundary vertex, we assume that each non-root vertex is associated with its 
parent in the tree. This allows us to find also u' in constant time. 

As for u", suppose we have precomputed, for each boundary vertex vj and 
each w & V \ {vj}, the successor of Vj on the path from Vj to w in shortest 
path tree T{vj). Depth-first searches in each shortest path tree allow us to 
make these precomputations in 0(n^/^) time and space. 

Now, since shortest paths are unique, the subpath of Pu from to m is a 
path in shortest path tree T{vi^) and u" is the successor of Vi^ on this path. 
With the above precomputations, we can thus find u" in constant time. 

Finally, to determine whether P„ is part of a clockwise walk of C, we do 
as follows. We first find the elementary faces adjacent to e in G. They can 
be obtained from dual tree T(f j) in constant time. We can also determine in 
constant time which of the two elementary faces belongs to int{G) since that 
elementary face is a child of the other in T{vi). Let fuv be the elementary 
face in the interior of G. We check if the edge directed from m to f is part of 
a clockwise or counter-clockwise walk of fuv Again, this takes constant time. 
If it is part of a clockwise walk of fuv then Pu is part of a clockwise walk of 
G (Figure [T2t^a)) and otherwise, Pu is part of a counter-clockwise walk of G 
(Figure [Hb)). 

This concludes the proof of Lemma [TOl 

Proof of Lemma [TT] 

Assume first that e is not an edge of Gi. Let Pi and P2 be the two shortest 
paths in G from v to the end vertices of e, respectively. Since e is not in 
Gi, it must belong to 6*2 • Hence, the intersection between G and Gi is the 
union of paths Q, where Q is a subpath of either Pi or P2 with both its end 
vertices in Vj. Each such path Q is a shortest path in Gi. It then follows 
from Lemma [6] that G does not cross any cycle of B[. 

Now, assume that e belongs to Gi and let /i, /2, and fj be defined as in 
the lemma. Let G' G B'l be given. We consider two cases: C intiC') and 
J C ext{G'). 
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Assume first that J C int{C'). Then R{C', Bi) is an ancestor of R{fj, Bi). 
Since vertex v of C belongs to Vj, part of C is contained in int{C'). 

It follows that if C does not cross C then e is contained in intiC'). The 
converse is also true. For if e is contained in intiC') then by Lemma [6|, both 
Pi and P2 are contained in int{C'), implying that C does not cross C. 

Thus, C crosses C if and only if e is not in int{C'), i.e., if and only 
if /i and /2 are both contained in ext{C'). The latter is equivalent to the 
condition that R{C',Bi) is an ancestor of neither R{fi,Bi) nor R{f2,Bi) in 
T{Bi). Hence, C crosses C if and only if the second condition of the lemma 
is satisfied. 

Now, assume that J C ext{C'). Then R{C',Bi) is not an ancestor of 
R{fj,Bi). Again, Lemma E] shows that C crosses C if and only if e is not 
in ext{C'), i.e., if and only if /i and /2 are both contained in int{C'). This 
holds if and only if R{C', Bi) is an ancestor of both Bi) and R{f2, Bi) 
in T{Bi). It follows that C crosses C if and only if the first condition of the 
lemma is satisfied. 

Proof of Lemma [12] 

Assume first that C ext{C) and let C be a descendant of C in T{Bi). 
We need to show that C is added to B. Since C G there is a pair of 
elementary faces /i and /2 in Gi which are separated by C and not by any 
other cycle in Bi. Let /i be contained in int(C') and let /2 be contained 
in ext{C'). Note that /2 is contained in int{C) since otherwise, C would 
separate /i and /2. 

Since JT" C ext{C) and since no cycle of B crosses C, all cycles oi B\ B[ 
belong to ext{C). Hence, no cycle of ^{Vj) U B[ U i32 \ {C} separates /i and 
/2. Since the set of cycles in the GMCB of G is a subset of ulVj) U B[ U B'^ 
by Lemma m and since C E B[, it follows that C is added to B. 

Now assume that JT" C int{C). Since no cycle of B crosses C, all cycles of 
B\B'i belong to int{C). A similar argument as the above then shows that 
all cycles of Bi belonging to ext{C) must be part of the GMCB of G. These 
cycles are exactly the those that are not descendants of C in T{Bi). 
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